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ANALYSIS  OF  HEAT-TOANSFER  EFFECTS  IN  ROCKET  NOZZLES 

OPERATING  WITH  VERY  HIGH-TEMPERATURE  HYDROGEN 

> 

■by  John  R.  Howell^  Mary  K.  Strite^  and  Harold  E.  Rehkel 
Lewis  Research  Center 


SUMMARY 

I  An  analytical  technique  suitable  for  the  solution  bf  complex  energy  trans¬ 
fer  problems  involving  coupled  radiant  and  convective  energy  transfer  is  de¬ 
veloped.  Solutions  for  the  coupled  axial  -wall  energy  flux  distribution  in  roc¬ 
ket  nozzles  using  hydrogen  as  a  propellant  are  presented.  Flow  rates  and  tem¬ 
peratures  studied  are  near  those  forecast  for  gaseous-core  nuclear -propulsion 
systems.  Parameters  varied  are  nozzle  shape ^  iiolet  propellant  temperature^ 
mean  reactor  cavity  temperature^  and  nozzle  wall  temperature  level.  The  effects 
of  variation  of  the  propellant  radiation  absorption  coefficient  with  pressure^ 
temperature^  and  wavelength  are  presented,  and  real  property  variations  are 
used  where  they  appear  to  be  significant.  Comparison  is  made  to  a  simplified, 
coupled  solution  using  a  modified  second-order  one-dimensional  diffusion  eq_ua- 

tion  for  the  radiative  transfer,  f  J 

. i  ^  ' 

At  the  temperature  levels  assumed,  radiative  transfer  may  account  for  a 
greater  portion  of  the  total  energy  transfer  over  important  portions  of  the 
nozzle,  and  its  effects  cannot,  therefore,  be  neglected. 

Extreme  energy  fluxes  (near  3X10®  Btu/(hr)(sq,  ft))  are  observed  for  cer¬ 
tain  cases,  and  this  implies  that  new  nozzle  cooling  techniq.ues  must  be  de¬ 
veloped. 
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With  the  advent  of  advanced  nuclear-propulsion  systems  such  as  the 
gaseous -core  configuration,  extremely  high  propellant  temperatirr.es  are  ex¬ 
pected.  Propellant  temperatures  encountered  in  the  nozzles  of  such  systems 
may  exceed  15,000^,  leading  to  the  belief  that  radiative  energy  transfer  in 
the  nozzle  may  be  of  significant  importance  in  comparison  with  convective 
heat  transfer.  The  purpose  of  this  report  is  the  investigation  of  the  relative 
importance  of  these  modes  of  energy  transfer  at  the  temperature,  levels  expected 
and  the  development  of  a  suitable  analytical  method  for  examing  them.  J 

Themal  design  of  a  nozzle  depends  on  knowledge  of  the  axial  heat-flux"^'' 


distribution  at  the  nozzle  wall.  The  flux  at  any  point  along  the  nozzle  is  com¬ 
posed  of  convective  and  radiative  contributions .  The  convective  heat  transfer 
is  affected  by  variations  in  propellant  enthalpy,  viscosity,  Prandtl  number, 
Mach  number,  specific  heat  ratio,  and  molecular  weight,  each  of  which  is,  in 
turn,  a  function  of  local  conditions.  Radiative  energy  transfer  is  dependent 
on  thfe  distribution  of  propellant  spectral  absorption  coefficient  for  radiation, 
which  is  dependent  on  the  local  static  temperature  and  pressure  distribution 
throughout  the  nozzle .  Both  radiative  and  convective  energy  transfer  are 
coupled  to  the  propellant  flow.  All  of  these  factors  are  also  geometry  depen¬ 
dent  . 


Robbins  (ref.  l)  and  Robbins,  Bachkin,  and  Medeiros  (ref.  2)  studied  the 
thermal  design  of  a  nozzle  in  conjunction  with  a  solid-core  nuclear  reactor. 
Both  papers  considered  radiant  transfer  between  surfaces  within  the  nozzle,  but 
the  propellant  was  assumed  to  be  transparent,  and  absorption  and  emission  of 
radiant  energy  in  the  propellant  were  therefore  neglected.  Coupling  between 
radiant  and  convective  heat  transfer  in  the  propellant  does  not  occur  in  this 
case.  G-rueber  (ref.  3)  studied  a  nozzle  cooled  solely  by  radiation  but  also 
assumed  the  propellant  transparent.  For  some  propellants,  notably  hydrogen, 
assumption  of  transparency  at  the  temperature  levels  considered  previously 
(below  6000°  R)  is  good;  however,  it  becomes  progressively  poorer  as  the  tem¬ 
perature  levels  increase . 

Convective  heat  transfer  in  rocket  nozzles  has  received  considerable  at¬ 
tention  because  it  is  the  dominant  heat -transfer  mode  in  the  systems  widely 
used  today.  Representative  papers  are  those  of  Bartz  (refs.  4  and  5),  Fortini 
and  Ehlers  (ref.  6),  Benser  and  Graham  (ref.  7),  and  Welsh  and  Witte  (ref.  8). 

Neumann  and  Bettinger  (ref.  9)  examined  the  major  accepted  methods  of  com¬ 
puting  the  gas-side  convective  heat-transfer  coefficients  in  solid-core 
nuclear-rocket  systems  and  demonstrated  the  effects  of  many  variables  on  the 
predicted  convective  heat  flux. 

Of  interest,  when  the  coupling  of  radiation  with  flow  is  considered,  are 
a  number  of  analyses  limited  chiefly  by  their  restriction  to  simple  geometries 
and  constant  properties.  Some  of  those  that  most  closely  approximate  the 
problems  studied  herein  are  as  follows:  Einstein  (refs.  10  and  11 )  obtained 
heat-transfer  rates  and  temperature  distributions  between  semi-infinite  plates 
and  in  circular  ducts  containing  flowing  gases.  He  considered  the  effects  of 
flow,  conduction,  and  radiation  in  the  gray  absorbing-emitting  medium,  and 
included  Internal  energy  generation  as  a  parameter.  Chen  (ref.  12),  using  an 
approximate  formulation,  was  able  to  study  radiation-scattering  effects  under 
similar  conditions  between  semi-infinite  plates. 

Kramer  (ref.  13)  examined  combined  convective  and  radiative  transfer  for 
use  in  an  analysis  of  transpiration  cooling,  but  he  considered  the  propellant 
to  be  a  "black"  gas  radiating  as  a  black  surface  directly  to  the  nozzle  wall 
and  did  not  include  coupling  between  radiation  and  convection.  No  considera¬ 
tion  was  given  to  the  high-temperature  reactor  chamber  or  its  contribution  of 
radiant  energy  to  the  nozzle. 
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Howell,  Strite,  and  Renkel  (ref.  14)  analyzed  a  specific  nozzle  shape  for 
the  coupled  convective  and  radiative  transfer.  The  first  paper  was  restricted 
to  constant  propellant  properties  with  a  cursory  examination  of  convective  ef¬ 
fects,  while  the  second  included  variations  in  propellant  properties. 

This  report  examines  in  detail  the  radiative  contribution  to  the  total 
axial  heat-flux  distribution  in  a  nozzle.  Property  variations  and  their  effect 
on  the  radiative  and  convective  heat  flux  and  radiation  from  the  upstream  reac¬ 
tor  chamber  are  taken  into  account . 


SYMBOLS 

Aj^  area  of  nozzle  wall  bounding  increment  i,  sq.  ft 
Ax  cross-sectional  area  of  nozzle,  sq.  ft 
B  convergence  control  constant 

Cg  energy  per  bundle,  Btu 

C-n  i  heat  capacity  of  propellant  evaluated  at  conditions  in  Increment 
i,  Btu/Clb)^^) 

D  local  nozzle  diameter,  ft 

e-^  Planck  spectral  distribution  of  emissive  power,  Btu/(sq.  ft)(hr)(ft) 
Pjj(r)  radial  emission  function  at  nozzle  entrance,  ft 
^  radiation  exchange  factor 

H  propellant  enthalpy,  Btu/ lb 

h^  nozzle  wall  heat-transfer  coefficient  in  increment  i,  lb/(hr)(sq  ft) 
i  axial  increment  index 

j  radial  increment  index 

k  dummy  axial  increment  index 

L  path  length  between  radiating  and  absorbing  element,  ft 

Z  actual  bundle  path  length,  ft 

M  local  propellant  Mach  number 

m  dummy  radial  increment  index 
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N  net  rate  of  energy  ‘bundles  entering  nozzle^  ‘bundles/hr 

n  ‘bundle  index  for  radiant  source 

static  pressure  in  increment  i,  lb/sq_  ft 
Pr  propellant  Prandtl  number 

p  bundle  index  for  flow  source 

Q  rate  of  radiant  energy  entering  nozzle  from  reactor  chamber,  Btu/hr 

R  number  chosen  at  random  from  set  of  numbers  evenly  distributed 

in  range  0  to  1 

R  local  radius  of  nozzle^  ft 

r  ralial  position  from  nozzle  axis  of  bundle  at  end  of  path^  ft 

rad  sq_uare  of  radial  bundle  position  at  end  of  path;  r^;  sq  ft 

S  rate  of  bundle  absorptions  in  nozzle  vail  increment;  bundles/hr 

S(i;j)  rate  of  bundle  absorptions  in  propellant  increment  (i;j); 

bundles/hr 

T  temperature;  ^ 

U(i^j)  rate  of  bundle  emissions  due  to  flow  sources  in  increment  (i;j), 
bundles/hr 

V  volume;  cu  ft 

W  total  propellant  mass-flov  rate;  Ib/hr 

Wj  propellant  mass-flov  rate  in  increment  j;  Ib/hr 

V  molecular  veight 

X  coordinate  perpendicular  to  nozzle  axis  (taken  through  bundle 

origin);  ft 

y  coordinate  perpendicilLar  to  x  and  Z;  ft 

z  axial  coordinate;  ft 

ai  angle  betveen  nozzle  vail  and  nozzle  axis  at  increment  i 

P  program  index:  0  for  radiant  source  segment;  1  for  flov  source 

segment 
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(x^  +  y2)  where  x  and  y  are  computed  at  position  of  interest^ 
sq_  ft 

propellant  ratio  of  specific  heats 

desired  percentage  temperature  change  between  iterations 
angle  measured  clockwise  from  x  around  z 
spectral  propellant  radiation  absorption  coefficient^  ft“^ 
mean  propellant  radiation  absorption  coefficient ^  ft"^ 
wavelength^  ft 

propellant  viscosity^  lb/(ft)(hr) 
radiation  path  length  in  volume  element  dV  ^  ft 
Stefan-Boltzmann  constant^  1.714X10  ^  Btu/(hr)(sq_  ft)(  R  ) 
Subscripts: 
ad  adiabatic 

a^e  absorbing  or  emitting  element 

i  centerline 

ch  chamber 

conv  convective 

d  truncated  at  decimal  point 

g  gas 

i  evaluated  in  increment  i 

in  inlet 

i^j  evaluated  in  increment  i^  3 

3  evaluated  in  increment  3 

max  maximum 

min  minimum 


r 

r 

5 

^  K 

K 

* 

A 

I 

o 
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n  nozzle 

nev  results  of  cnxrent  iteration 
o  evaluated  at  origin 

rad  radiated  portion 

ref  evaluated  at  reference  condition 
s  static  condition 

St  stagnation 

t  evaluated  at  nozzle  throat 

w  vail 

A  wavelength  dependent 


WHOD  OF  MALYSIS 

The  analytical  method  consists  of  writing  an  energy  balance  on  each  volume 
element  of  propellant  within  the  nozzle  shown  in  figure  1  and  solving  the  re- 


Solid  lines  bound  volume 
where  energy  fluxes  and 
temperatures  are 
calculated^ 


Figure  1.  -  Model  for  nozzle  heat  transfer. 
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sultant  set  of  norainear  second-order  integro-differential  equations.  A  typical 
equation  from  this  set  and  a  discussion  of  the  detailed  technique  used  in  the 
solution  are  given  in  appendix  A. 

A  complete  flow  scheme  showing  the  application  of  these  equations  to  a 
digital  computer  program  is  given  in  appendix  B.  A  discussion  of  computer 
running  time  and  program  convergence  is  included. 

In  hrief^  the  method  of  solution  involves  evaluation  of  the  integrals  im¬ 
plicit  in  the  equations  hy  a  Monte  Carlo  technique.  A  new  axial  total  temper¬ 
ature  distribution  is  then  found  by  solving  the  set  of  simultaneous  equations 
by  the  Newton-Raphson  method  as  presented  in  reference  15.  The  axial  static 
temperature  distribution  is  then  computed  and  used  to  find  local  properties. 

This  prodedure  is  repeated  until  convergence  of  the  temperature  distribution. 
Axial  heat-flux  distributions  are  then  computed  on.  the  basis  of  this  temperature 
profile . 


Assumptions 

In  order  to  make  the  analysis  tractable^  the  following  assimiptions  are 

made: 


(1)  The  nozzle  walls  are  perfect  absorbers  for  all  incident  radiation. 

(2)  No  energy  is  added  to  the  propellant  by  radiant  emission  from  the 
nozzle  walls. 

(3)  The  radiation  absorption  coefficient  in  the  propellant  does  not  vary 
along  the  path  between  the  points  of  energy  emission  and  absorption^  and  its 
value  is  based  on  local  conditions  at  the  point  of  emission. 

(4)  A  one -dimensional  isentropic  analysis  is  adequate  to  describe  the  Mach 
number  and  the  static  pressure  distributions  in  the  nozzle. 

(5)  Equilibrium  physical  properties  can  be  used  at  local  conditions. 

(6)  Gamma  heating  of  the  propellant  and  the  nozzle  walls  can  be  ignored. 

(7)  Complete  radial  mixing  occurs  within  the  nozzle^  so  that  there  are  no 
radial  temperature  gradients . 

(8)  No  net  flow  crosses  one -dimensional  streamlines  within  the  nozzle. 
These  assumptions  and  their  validity  are  discussed  in  appendix  A. 


Features  of  - Method 

Incorporation  of  the  Monte  Carlo  technique  to  evaluate  the  radiation  terms 
allows  freedom  from  the  common  transparent  or  constant  absorption  coefficient 
propellant  assumptions  and  further  allows  a  reasonably  straightforward  analysis 
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of  the  coupling  between  radiation  and  convection.  The  complete  solution  tech¬ 
nique  allows  consideration  of  most  variations  of  physical  properties  with  tem¬ 
perature  and  pressure  where  these  variations  appear  important.  As  finally  pro¬ 
grammed  for  the  digital  computer^  these  effects  plus  those  of  nozzle  geometry^ 
mass  flow  rate  and  mass  flow  rate  distribution^  reactor  cavity  radiation^  and 
nozzle  wall  temperature  distribution  can  be  investigated. 


Diffusion  Solution 

A  much  more  straightforward^  but  also  more  restricted^  solution  to  the 
nozzle  heat-transfer  problem  can  be  obtained  by  modifying  the  second-order  dif¬ 
fusion  solution  as  presented  by  Deissler  (ref.  16)  for  the  radiative  portion  of 
the  problem  and  combining  this  with  the  pipe  flow  equation  for  convection.  The 
assumptions  involved  and  the  method  of  modification  used  to  obtain  these  equa¬ 
tions  are  given  in  appendix  C. 

It  is  possible  to  obtain  a  solution  by  hand  with  this  method^  whereas  a 
high-speed  digital  computer  is  mandatory  for  the  Monte  Carlo  solution.  Neither 
the  effects  of  thermal  radiation  from  the  nuclear -reactor  cavity  nor  the  spec¬ 
tral  effects^  however^  are  conveniently  studied  by  the  one -dimensional  radiation 
diffusion  approach  used  herein. 

Determination  of  the  heat  transfer  is  obtained  by  iterative  solution  of 
the  equation 
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starting  with  volunie  element  1=1.  Solution  for  this  element  gives  an  initial 
guess  for  element  i  =  2,  and  this  procedure  is  continued  along  the  nozzle 
length. 

The  axial  heat -flux  distribution  is  then  obtained  from 


WC 

~ 


(2) 


DESCRIPTION  OF  PROBLEM  ANALYZED 

The  problem  solved  by  the  methods  presented  is  the  heat -transfer  rate  as  a 
function  of  axial  position  for  a  conical  nozzle.  The  propellant  is  assumed  to 
be  hydrogen.  The  effects  of  nozzle  shape,  radiation  absorption  coefficient, 
inlet  propellant  temperature,  apparent  reactor  cavity  temperature,  variable 
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Length,  ft 

Figure  2.  -  Dimensions  of  nozzles. 


TABLE  I.  -  MASS  FLOW  RATES  OF  HTDROGEN 
FOR  VARIOUS  NOZZLE  SHAPES 


nozzle  wall  temperature,  and  the  coupling  of 
radiation  and  convection  will  be  demon¬ 
strated. 

The  effects  of  a  number  of  parameters 
on  the  axial  distribution  of  energy  flux  to 
the  walls  of  a  nozzle  are  presented.  For 
each  nozzle  shape  examined  (fig.  2),  the 
flow  rate  of  propellant  was  taken  as  that 
producing  sonic  flow  at  the  throat  on  the 
basis  of  the  one -dimensional  isentropic  flow 
eq_uations .  These  flow  rates  are  shown  in 
table  I.  The  nozzle  shapes  and  sizes  were 
chosen  as  representative  of  those  anticipated 
for  gaseous -core  nuclear  systems. 


[Entering  propellant  static  pressure^ 
100  atm.] 


Nozzle 

Inlet  temperature;  ^ 

15; 000 

15; 000 

Mass  flow  rate  of  hydrogen;  Ih/ sec 

I 

940 

— 

II 

350 

320 

III 

153 

Propellant  enthalpy,  viscosity,  heat 
capacity,  ratio  of  specific  heat,  and  radia¬ 
tion  absorption  coefficient  were  allowed  to 
vary  in  the  manner  noted  in  appendix  C  for 
all  the  results  presented.  On  the  basis  of 
the  one -dimensional  flow  assumption,  the 
propellant  is  divided  into  volume  elements; 
typical  ones  are  shown  in  figure  1. 

Pressure  and  Mach  number  distributions 
calculated  from  standard  isentropic  flow 
eq^uations  are  assimied  sufficiently  accurate 
and  are  shown  in  figures  3  and  4,  respec¬ 
tively,  for  the  nozzle  shapes  studied.  Ex¬ 
cept  for  the  figure  showing  the  effect  of 
wall  temperature,  all  calculations  are  based 
on  a  constant  nozzle  wall  temperature  of 
4000°  R. 


RESULTS  AND  DISCUSSIOE 
G-eneral  Remarks 

A  cursory  examination  of  figures  5  to  11 
shows  that  the  shape  of  the  axial  energy  flux 
distribution  in  the  nozzle  is  similar  for  all 
cases  where  radiation  and  convection  are  con¬ 
sidered  simultaneously.  A  large  peak  in 
total  flux,  reaching  values  near  2X10®  Btu/ 
(hr)(sq_  ft),  occurs  at  the  nozzle  throat  be¬ 
cause  of  the  increased  mass -flow  rate  per 
unit  area. 

The  total  flux  drops  rapidly  downstream 
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Figure  4.  -  Nozzle  Mach  number  distributions. 


0  1  2  3  4  5  6  7  8  9 

Axial  position,  ft 


Figure  5.  -  Effect  of  radiation  absorption  coefficient  on  total 
radiant  energy  flux.  Nozzle,  11;  chamber  temperature, 
13,000®  R;  inlet  temperature,  13,000®  R;  wall  tempera¬ 
ture,  4000®  R. 


of  the  throat  because  of  a  number 
of  interrelated  effects.  The  total 
propellant  temperature  has  de¬ 
creased  slightly  because  of  up¬ 
stream  heat  losses^  and  the  propel¬ 
lant  static  temperature  is  also 
lowered  rapidly  at  points  further 
downstream  because  of  the  increas¬ 
ing  propellant  Mach  number.  Both 
of  these  effects  tend  to  lower  the 
convective  heat  transfer. 

The  radiative  transfer^  which 
can  be  quite  significant  upstream 
of  the  nozzle  throaty  also  de¬ 
creases  rapidly  downstream  of  the 
throat.  This  is  because  the  fourth 
power  of  the  static  temperature  is 
decreasing  rapidly^  and  the  hydro¬ 
gen  propellant  becomes  nearly 
transparent  at  the  lower  static 
conditions  downstream  of  the 
throat.  Both  factors  tend  to  radi¬ 
cally  decrease  the  radiative  trans¬ 
fer.  Radiation  contributes  from 
under  5  to  nearly  90  percent  of  the 
total  flux  depending  on  the  case 
under  consideration. 


Effect  of  Variation  of  Radiation 

Absorption  Coefficient 

Figure  5  demonstrates  the  ef¬ 
fects  on  the  radiative  energy  fl-ux 


of  certain  assumptions  about  the 

variation  of  propellant  absorption  coefficient  with  local  conditions.  While  the 


constant  mean  absorption  coefficient  chosen  predicts  a  greater  radiative  trans¬ 
fer^  the  other  assumptions  lead  to  results  that  are  in  fair  agreement  except 
in  the  divergent  section  where  radiation  effects  are  small.  This  was  the  basis 


for  computing  most  of  the  results  by  using  only  the  simplest  variation^  that 
is^  with  static  pressure.  The  assumption  of  variation  with  both  static  temper¬ 
ature  and  pressure  predicted  the  lowest  radiative,  transfer  rate  as  expected 
from  the  property  variations  for  this  system  shown  in  figure  4.  Since  each 
variation  examined  lowered  the  computed  flux  from  that  found  for  a  constant 
coefficient,  consideration  of  variation  with  all  three  (static  temperature, 
pressure,  and  wavelength)  simultaneously  would  be  expected  to  predict  even 
smaller  radiant  fluxes. 
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Nozzle 

1  .  . 

1 

;in 

1  7 

'hroa' 

\ 

\ 

J 

i\. 

■ 

\I 

1  11^ 

1  2  3  4  5  6  7  8  9 

Axial  position,  ft 

Figure  6.  -  Effect  of  nozzle  shape  on  independent 
radiant  flux  to  nozzle.  Inlet  temperature, 

13, 000®  R;  reactor  temperature,  13, 000®  R; 
wall  temperature,  4000®  R;  i<  =  F(P>y=i3  ooqO  r* 


0  1  2 


3  4  5  6 

Axial  position,  ft 


Figure  7.  -  Independent  convective  flux  to  nozzle. 
Inlet  temperature,  13, 000®  R;  wall  temperature, 
4000®  R. 


Effect  of  Nozzle  Shape  on  Flux  Distribution 

In  figure  6^  the  radiant  flux  is  shown  with  the  nozzle  shape  as  a  param¬ 
eter.  Here  the  absorption  coefficient  is  assumed  to  be  dependent  only  on  local 
static  pressure.  Since  it  has  the  largest  throat  diameter  and;  therefore^  the 
smallest  convergence  angle;  nozzle  I  has  the  smallest  static  pressure  and  tem¬ 
perature  drop.  This  in  turn  causes  it  to  have  the  highest  radiative  flux  past 
the  throat  because  the  absorption  coefficient  remains  comparatively  high;  as 
does  the  local  static  temperature.  In  addition;  geometric  shielding  is  less. 

On  the  other  hand;  the  higher  Mach  number  and  lower  local  static  temperature  and 
pressure  at  a  given  axial  location  upstream  of  the  throat  in  nozzle  I  results  in 
a  comparatively  lower  local  radiative  flux  in  this  part  of  the  nozzle. 

Figure  7  indicates  the  effect  of  nozzle  shape  on  the  convective  heat  trans¬ 
fer.  Nozzle  III  has  the  smallest  throat  diameter  and;  thuS;  the  greatest  con¬ 
vergence  angle  and  the  smallest  flow  rate  of  propellant.  This  causes  the  con¬ 
vective  flux  in  the  convergent  and  divergent  portions  of  nozzle  III  to  be  smal¬ 
ler  than  that  in  the  other  nozzles.  At  the  throat;  however;  nozzle  III  exhib¬ 
its  the  largest  heat  flux;  this  is  due  to  the  effect  of  the  large  mass-flow  rate 
per  unit  area  on  the  heat-transfer  coefficient. 


Effect  of  Considering  Total  Heat  Fliix  as  Additive  or  Coupled 

The  radiative  fluxes  of  figure  6  and  the  convective  fluxes  of  figure  7 
were  computed  independently.  However;  because  each  heat-transfer  mechanism 
tends  to  lower  the  total  propellant  temperature;  which  in  turn  affects  the 
static  temperature  and  the  properties;  there  is  some  q.uestion  as  to  whether 
simply  adding  the  results  will  accurately  predict  the  total  heat  transfer. 
The  computer  program  is  written  to  consider  the  two  modes  of  heat  transfer 


Energy  flux,  Btu/IhrMsq  ft) 


Axial  position,  ft 

(a)  Comparison  of  coupled  and  independent  solu-  (b)  Comparison  of  coupled  and  diffusion  solutions, 
tions.  Chamber  temoerature,  13, 000°  R;  inlet 
temperature,  13, 000°  R;  wall  temperature, 

4000°  R. 


Figure  8.  -  Comparison  of  total  energy  flux,  k  -  F(P)^„i3  qqqO 


Figure  9.  -  Effect  of  inlet  propellant  temperature 
on  radiative  and  total  coupled  energy  flux. 
Nozzle,  II;  -  T|p;  wall  temperature, 

4000°  R;  K  •  F(P)|-.^. 


as  coupled^  and  this  is  done  in  the  remainder 
of  the  figures. 

Figure  8(a)  shows  the  differences  in  the 
total  energy  fluxes  obtained  hy  adding  the 
radiative  and  convective  fluxes  and  those 
resulting  from  calculating  the  coupled  ef¬ 
fect.  The  additive  solutions  are  seen  to 
slightly  overpredict  the  total  flirx,  because 
each  is  based  on  a  propellant  temperature 
distribution  that  is  somewhat  too  high. 

Since  the  energy  loss  in  the  propellant 
stream  is  relatively  small  at  the  large  flow 
rates  present  in  these  nozzles,  the  total 
temperature  of  the  stream  remains  almost  un¬ 
changed  no  matter  which  method  is  used  in 
calculating  the  heat  transfer.  The  use  of 
coupled  calculations  for  the  specific  cases 
presented  herein  is,  therefore,  not  neces¬ 
sary,  since  an  additive  method  produces  ac¬ 
ceptable  accuracy. 


Compartson  With  Solution  by  Radiation  Diffusion  Method 

Comparison  of  the  more  complete  solution  is  made  with  the  second-order  dif¬ 
fusion  solution  with  Jump  boundary  conditions  of  reference  17  and  modified  as 
described  in  appendix  C.-  The  diffusion  results  (fig.  8(b))  presented  are  based 


on  the  assumption  that  the  radiation  is  proportional  to  the  fourth  power  of  the 
total  temperature  of  the  propellant  rather  than  the  static  temperature^  and  use 
a  pressure  dependent  Rosseland  mean  radiation  absorption  coefficient  (see  appen¬ 
dix  C)  taken  from  reference  17.  Propellant  heat  capacity  was  allowed  to  vary 
with  pressure,  but  the  propellant  ratio  of  specific  heats  was  taken  as  constant 
and  equal  to  1.0.  This  caused  neglect  of  the  effects  of  static  temperatiire 
variation  due  to  acceleration  of  the  propellant. 

Convection  was  computed  as  for  the  more  complete  solution.  Comparison  of 
the  solutions  for  similar  cases  is  shown  in  figure  8(b).  As  expected,  the  more 
complete  solution  gives  more  accurate  values  of  energy  flux  upstream  of  the 
throat  because  the  energy  entering  from  the  upstream  reactor  chamber  is  con- 
Downstream  of  the  throat,  the  diffusion  solution  overpredicts  the 
flimc  because  variations  of  the  static  temperature  were  neglected.  Aside  from 
these  effects,  agreement  is  surprisingly  good.  For  certain  cases,  especially 
where  large  effects  of  the  reactor  chamber  thermal  radiation  are  present,  less 
agreement  would  be  expected. 


Effect  of  Inlet  Propellant  Temperature 

Figure  9  demonstrates  the  effect  of  two  propellant  temperatures  at  the  noz¬ 
zle  inlet  on  the  total  coupled  energy  flux.  As  expected,  the  flux  is  much 
greater  at  the  higher  propellant  temperature.  Two  factors  mainly  account  for 
this.  They  are  the  increased  temperature  level  of  the  propellant  which  in¬ 
creases  the  radiant  flux,  and  the  larger  enthalpy  difference  between  the  propel¬ 
lant  and  the  wall,  which  increases  the  convective  flux.  Another  factor  is  the 
property  variations  that  occur  with  temperature;  in  most  cases,  the  property 
variations  lend  an  additional  increase  to  the  energy  flux  at  higher  tempera¬ 
tures,  although  variations  in  the  heat  capacity  and  ratio  of  specific  heats  are 
so  complex  as  to  forestall  a  definite  statement  in  their  cases. 


Effect  of  Reactor  Chamber  Temperature  Level 

The  upstream  reactor  chamber  may  have  a  large  effect  on  the  nozzle  heat 
transfer  through  its  thermal  radiation  contribution.  Figure  10(a)  demonstrates 
the  magnitude  of  possible  variation,  based  on  a  number  of  mean  reactor  chamber 
temperatures.  The  radiative  contribution  to  the  energy  flux  becomes  so  large 
for  high  reactor  chamber  temperatures  that  the  peak  flux  occurs  at  the  nozzle 
entrance  rather  than  the  throat,  and  a  severe  cooling  problem  appears  at  this 
rather  unexpected  area. 

Since  the  reactor  radiation  acts  as  an  energy  source,  the  propellant  tem¬ 
perature  is  expected  to  change  locally  for  different  reactor  temperatures. 
Figure  10(b)  indicates  that  the  stagnation  temperature  of  propellant  within  the 
nozzle  may  increase  because  of  this  effect. 

Because  of  the  possible  wide  variation  between  the  temperatures  of  the 
major  emitting  (reactor  chamber)  and  absorbing  (nozzle  propellant)  portions  of 
the  system,  spectral  effects  could  become  important.  Assuming  the  propellant 
properties  to  be  independent  of  wavelength  leads  to  the  differences  in  calcu- 
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Energy  flux,  Btu/{hr)(sq  ft) 


temperature 


lated  propellant  temperature  shown  in 
figure  10(1)).  As  the  wavelength  effect  in¬ 
creases  the  propellant  temperature,  more 
energy  must  he  absorbed  by  the  propellant 
than  for  the  wavelength  independent  case. 
This  implies  that  less  energy  will  reach 
the  nozzle  walls  for  the  spectral  case,  and 
this  conclusion  is  borne  out  by  the  resiuLts 
presented  in  figure  10(c). 


Effect  of  Reactor  Chamber  Temperature 
Distribution 


1  2 


3  4  5  6 

Axial  position,  ft 


Figixre  10(d)  shows  the  effect  of  nonuni 
formlty  in  the  radial  temperature  distri¬ 
bution  in  the  reactor  cavity.  The  solid 
lines  show  the  energy  flux  profiles  for  a 
constant  reactor  cavity  temperature  of 
20,000°  R,  while  the  dotted  lines  are  the 
profiles  for  a  parabolic  temperature  dis¬ 
tribution  in  the  reactor  cavity  that  radiates  the  same  total  energy  as  in  the 
constant  temperature  case.  This  distribution  had  a  peak  temperature  at  the 
nozzle  centerline  of  24,000°  R  and  a  minimum  at  the  nozzle  wall  of  15,000°  R. 
A  significant  reduction  in  the  peak  radiative  fliuc  at  the  nozzle  wall  is  evi¬ 
dent,  indicating  that  the  distribution  of  entering  energy  is  important. 


Figure  11.  -  Effect  of  nozzle  surface  temperature  on 
total  coupled  energy  flux.  Nozzle,  II;  chamber  tem¬ 
perature,  13,000®^  inlet  temperature,  13, 000°  R; 

^  ■''<'’>1=13, 000®  R- 


Effect  of  Nozzle  Surface  Temperature 

The  nozzle  surface  temperature  is  a  relatively  weak  parameter  (fig.  ll). 

A  large  (4000°  R)  change  in  the  surface  temperature  produces  only  a  relatively 
small  variation  in  energy  flux. 


CONCLUDING  REMARKS 


An  analytical  method  for  determining  the  axial  distribution  of  coupled 
radiative  and  convective  energy  flux  on  a  rocket  nozzle  wall  is  introduced. 

The  method  consists  of  using  a  Monte  Carlo  technique  for  finding  the  radiative 
transfer  in  conjunction  with  a  standard  matrix  solution  routine  for  the  set  of 
governing  equations.  Speed  of  convergence  and  attainable  accuracy  in  solving 
a  nozz^  heat-transfer  problem  is  found  to  be 

1  The  results  of  the  analysis  for  cases  representative  of  anticipated  con¬ 
ditions  in  a  nozzle  used  with  a  gaseous-core  nuclear-propulsion  system  are  pre¬ 
sented.  The  magnitude  of  the  flvixes  computed  reach  extremely  high  levels  and, 
for  some  cases,  is  near  3X10^  Btu/(hr)(sq_  ft)_jJ\This  value  is  more  than  thirty 
times  higher  than  the  peak  values  computed  by  B^ser  and  Graham  (ref.  7)  for  a 
representative  large  nozzle  in  a  chemical -rocket  using  tQrdrogen-oxygen  propel- 


These  high  fluxes  may  occur  near  the  nozzle  entrance  and  are,  therefore, 
effective  over  a  relatively  large  area  of  the  nozzle  wall,  in  contrast  to 
chemical  rockets  where  the  peak  flux  occurs  only  over  a  relatively  small  area 
near  the  throat.  This  implies  extreme  heat  loads  to  any  nozzle  coolant  system. 

1  The  effect  on  the  flux  of  increased  wall  temperature  is  so  small  that  noz¬ 
zle  materials  with  higher  temperature  capabilities  will  be  of  little  value  in 
solving  the  cooling  problem.  Some  results  indicate  that  consideration  of  ra¬ 
dial  temperature  profiles  may  reduce  computed  peak  wall  fluxes.  1, 

\ 

'  The  flux  at  the  nozzle  wall  due  to  either  convection  or  radiation  is  of 
such  magnitude  as  to  demand  cooling  methods  involving  film  cooling,  transpira¬ 
tion,  houndary -layer  seeding,  or  possibly  elimination  of  a  physical  nozzle  and 
use  of  electromagnetic  effects^  (  \  TM' 
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APPEKDIX  A 


ANALYSIS 

Method  of  Solution 

For  a  general  volume  element  dV  with  axial  position  i;hdex  i  and  radial 
position  index  hy  using  the  assumptions  listed  in  the  n^ction  Assumptions, 
an  energy  halance  gives 

(Energy  radiated  into  \ 

from  other  propellant  elements)  +  (Flow  energy  entering  dVi^j)  = 
and  the  reactor  cavity  / 

(Energy  radiated  from  dVj_^j)  +  (Flow  energy  leaving  dV^^j) 

I 

+  (Energy  lost  to  nozzle  wall  from  (Al) 

This  can  he  written  mathematically  hy  substituting  the  following  relations: 
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The  index  of  refraction  of  the  propellant  is  taken  as  constant  and  equal  to  1.0 
and  therefore  does  not  appear  in  the  radiation  terms . 

Solution  of  equation  (A1)  for  the  total  temperature  of  the  volimie  element 
is  desired^  hecause^  knowing  this^  the  energy  flux  to  the  wall  can  he  found. 
Solution  of  this  equation^  however^  demands  a  complete  knowledge  of  the  propel¬ 
lant  temperature  distribution  for  evaluation  of  the  integral  terms.  Because  of 
this  and  the  temperature  dependence  of  many  of  the  physical  properties^  an 
iterative  procedure  must  he  used.  A  total  temperature  distribution  T(i)  is 
assimied  and  is  used  to  evaluate  the  static  temperature  distribution  that  is 
used^  in  turn^  to  evaluate  the  temperature  dependent  properties.  The  set  of 
equations  for  all  (i,j)  is  then  solved  simultaneously  by  the  Newton-Raphson 
method  (ref.  15)  to  obtain  a  new  total  temperatirre  distribution.  This  distri¬ 
bution  is  used  as  the  new  guess  to  determine  a  static  temperature  distribution 
for  evaluating  properties^  and  the  procedure  is  continued  until  the  assimied 
and  computed  temperature  distributions  agree  within  prescribed  limits. 

Solution  of  the  equations  is  greatly  complicated  by  the  presence  of  termed 
in  equation  (Al) .  Some  further  simplifications  must  be  made  to  make  an  evalu¬ 
ation  of  the  integral  terms  tractable.  One  desirable  simplification^  if  it  is 
justifiable  in  the  case  under  discussion^  is  to  consider  K  constant  along  the 
path  L  between  the  radiating  and  absorbing  element  even  though  it  remains 
dependent  on  the  static  temperature  and  pressure  of  the  radiating  element. 

Under  this  restriction  the  exponential  terms  inside  the  integrals  reduce  to 
exp  (-/CgL)  where  is  a  function  of  the  local  conditions  in  the  emitting 

element . 

Term  thus  modified^  can  be  evaluated  by  a  Monte  Carlo  technique.  This 
technique  consists  of  dividing  the  energy  emitted  by  each  radiating  volume  or 
area  into  finite  bundles^  and  following  each  bundle  through  its  lifetime  of 
absorptions  and  emissions  until  final  loss  from  the  enclosed  nozzle  volume. 

This  loss  can  occur  by  means  of  the  bundle  passing  through  the  entrance  or 
exit^  by  absorption  at  the  nozzle  wall  and^  in  some  cases^  in  propellant  flow 
sinks.  Derivation  of  the  functional  relations  governing  the  events  along  the 
paths  of  these  bundles  is  given  in  reference  18  with  a  discussion  of  applica¬ 
tion  to  radiant  energy  transfer  problems. 


Discussion  of  Assumptions 

Obviously^  limitations  are  placed  on  the  accuracy  of  the  analysis  by  the 
assumptions  made.  The  effect  of  many  of  the  assumptions  depends  on  the  speci¬ 
fic  nozzle  and  propellant  studied. 

Convective  and  radiative  energy  is  assumed  to  be  removed  from  each  radial 
element  in  proportion  to  the  radial  distribution  of  flow.  This  corresponds  to 
the  case  of  complete  radial  mixing  when  a  one -dimensional  flow  is  assumed. 

The  total  temperature  is  computed  solely  on  the  basis  of  heat -transfer  ef¬ 
fects;  the  lowering  of  total  temperature  due  to  friction  is  neglected.  In  the 
divergent  portion  of  the  nozzle^  therefore^  the  radiation  and  convection  from 
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the  propellant  will  be  somewhat  overpredicted. 

A  constant  radial  temperature  is  assumed  on  the  basis  of  the  large  flow 
rates  and  mixing  anticipated  in  such  systems.  Little  Information  is  available 
on  radial  mixing  for  nozzle  geometries,  and  the  validity  of  this  assumption 
cannot  be  assessed. 

A  further  assumption  is  that  the  propellant  radiation  absorption  coeffici¬ 
ent  is  a  function  of  at  most  two  variables:  either  (l)  local  propellant  static 
temperature  and  pressure  or  (2)  wavelength  and  local  propellant  static  pressure. 
It  is,  in  principle,  possible  to  include  the  variation  of  absorption  coefficient 
with  wavelength  together  with  the  effects  of  temperature  and  pressure.  However, 
computer  running  time  would  be  significantly  Increased  by  the  addition  of ^an¬ 
other  variable,  so  the  program  was  limited  to  variations  orly  in  the  combina¬ 
tions  noted.  For  any  propellant,  the  assumption  of  the  absorption  coefficient 
being  invariant  between  the  point  of  emission  and  the  point  of  absorption  be¬ 
comes  valid  if  any  of  the  following  criteria  are  met:  (l)  Conditions  through¬ 
out  the  nozzle  vary  only  slightly.  (2)  Distances  between  regions  at  different 
conditions  are  large  compared  to  the  radiation  mean  free  path  (the  reciprocal 
of  the  absorption  coefficient).  (3)  The  absorption  coefficient  is  a  weak  func¬ 
tion  of  local  conditions.  (4)  The  geometry  is  such  that  points  under  widely 
different  conditions  cannot  see  one  another.  Under  these  conditions,  either 
the  mean  absorption  coefficient  will  be  nearly  constant  throughout  the  nozzle 
or  radiant  energy  will  be  blocked  or  greatly  attenuated  before  reaching  a  re¬ 
gion  with  a  different  absorption  coefficient. 

Under  conditions  (l)  and  (3),  neglect  of  variation  of  the  absorption  coef¬ 
ficient  with  wavelength,  static  pressure,  or  static  temperature  also  becomes 
quite  a  good  assumption.  Variations  from  these  restrictions,  of  course,  add 
increasing  uncertainty  to  the  magnitude  of  the  errors  involved.  For  hydrogen 
propellant  in  a  nozzle  studied  in  the  example  herein,  these  assumptions  are  met 
quite  well.  This  is  because  the  major  part  of  the  radiative  transfer  occurs  in 
the  portion  of  the  nozzle  upstream  of  the  throat  where  static  pressure  and  tem¬ 
perature  variations  are  small  due  to  the  relatively  small  changes  in  Mach  num¬ 
ber.  This,  in  turn,  causes  the  local  property  variations  in  this  region  to  be 
small.  The  downstream  portion  of  the  nozzle,  where  variations  in^propell^t 
absorption  coefficient  can  be  significant  due  mainly  to  large  variations  in 
local  static  conditions,  contributes  little  radiant  energy  to  the  upstream  sec¬ 
tion  because  of  the  blocking  effect  of  the  convergent -divergent  geometry. 

Since  the  inlet  propellant  temperatures  studied  are  above  10,000°  R  and 
the  nozzle  wall  temperature  is  below  6000°  R,  the  radiant  energy  contributed  to 
the  system  by  the  wall  is  negligible  in  comparison  with  the  energy  entering 
from  the  reactor  cavity  and  that  originating  in  the  propellant.  For  this  rea¬ 
son,  no  energy  is  considered  to  be  originating  at  the  nozzle  wall  for  the  pro¬ 
pellant  heat  balance  even  though  an  emission  term  is  included  in  the  energy 
balance  used  to  determine  the  net  wall  heat  flux. 

To  account  for  energy  radiated  into  the  nozzle  from  the  reactor  cavity,^ 
the  assumption  is  made  that,  in  the  plane  normal  to  the  nozzle  axis,  the  radi¬ 
ant  energy  enters  diffusely,  that  is,  follows  the  cosine  law.  The  amount  of 
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radiation  entering  the  nozzle  from  the  reactor  chamber  must  be  specified.  For 
a  spectral  solution^  the  spectral  energy  distribution  is  taken  to  be  that  of  a 
black  body  at  the  mean  reactor  chamber  temperature. 

The  exhaust  plume  of  the  nozzle  is  assumed  to  be  transparent  to  radiation 
on  the  basis  of  the  mean  absorption  coefficient  data  for  hydrogen  shown  in 
figure  4  and^  therefore^  neither  absorbs  nor  emits  radiation.  Because  nuclear 
propulsion  will  in  general  be  used  only  in  space^  no  radiation  from  atmospheric 
interactions  is  considered. 

For  computing  the  value  of  the  convective  heat-transfer  coefficient^  the 
method  presented  by  Benser  and  Graham  (ref.  7)  was  chosen  as  being  relatively 
simple  to ^ use  while  being  reasonably  accurate.  The  definition  of  reasonably 
accurate  is  somewhat  nebulous  in  this  case  since  no  experimental  data  is  avail¬ 
able  under  the  conditions  imposed^  however^  at  lower  temperatures^  experimental 
comparison  seems  adequate  (ref.  6)  for  a  similar  analysis. 

Use  of  a  possibly  more  accurate  but  definitely  much  more  complex  boundary- 
layer  analysis  for  the  convective  heat  transfer  does  not  seem  justified  for 
this  work.  Certainly  the  pipe-flow  equation  as  modified  should  accurately  re¬ 
flect  the  effects  of  parametric  variations  in  the  significant  variables . 

Benser  and  Graham  (ref.  7)  suggest  the  use  of  an  enthalpy  difference  heat- 
transfer  coefficient  given  by  the  following  equation: 


0.026 


h  = 


W) 


^0.2  „  0.667 
D  Prref 


lb/ (hr) (sq  ft) 


with  the  properties  evaluated  at  a  reference  enthalpy  given  by 
Hpef  =  Hg  +  0.5  (Hj^  -  Hg)  +  0.22  Pr^^^  (H  -  Hg) 
The  convective  heat  flux  is  then  given  by 


(A2) 


)  -  -  V 

conv 

where 

Had  =  Hg  +  Pr^/|  (H  -  Hg) 

The  coefficient^  0.026,  in  equation  (A2)  is  an  arbitrary  value  that  baR 
been  shown  experimentally  to  vary  with  the  axial  position  in  the  nozzle.  Such 
variation  must  be  determined  experimentally  for  a  specific  nozzle  and  specific 
conditions  and  was  therefore  taken  as  constant  herein.  The  value  0.026  is  sug' 
gested  in  reference  7. 
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APPEEDIX  B 


COMPUTER  PROGRAM 

The  SHnplified  flow  chart  given  hy  figure  12(a)  shows  in  a  general  form  the 
computer  approach  followed  in  the  determination  of  the  nozzle  heat -flux  distri¬ 
bution.  The  major  sections  of  the  program  are  expanded  to  give  a  complete  flow 
chart  in  the  subsequent  parts  of  figure  12. 

Running  time  for  the  program  varied  widely .  As  finally  modified  and  pro¬ 
gramed  using  the  Fortran  IV  compiler  on  an  IBM  7094  digital  computer ^  complete 
running  times  for  a  single  case  averaged  much  less  than  10  minutes.  This 
varied  depending  on  the  astuteness  of  the  original  guess  of  propellant  tempera¬ 
ture  distribution,  the  magnitude  of  property  variations  in  the  static  pressure 
and  temperature  range  under  consideration,  the  accuracy  of  solution  required, 
and  a  number  of  other  interrelated  factors.  The  Monte  Carlo  segment  of  the 
program  was  the  largest  time  consumer,  and  the  time  used  was  especially  depend¬ 
ent  on  the  number  of  energy  bundles  followed  and  the  propellant  opacity. 


Accuracy  of  the  solutions  is  almost  impossible  to  bracket  analytically  in 
an  iterative  Monte  Carlo  solution  tied  to  other  finite  difference  program  seg-^ 
ments.  A  discussion  of  this  is  given  in  reference  19.  Lacking  such  an  analyti¬ 
cal  estimate,  however,  it  is  nevertheless  possible  to  deduce  the  error  by  more 
prosaic  means. 

For  these  solutions,  accuracy  could  be  determined  by  comparison  to  known 
limiting  solutions,  such  as  for  pure  convection  or  pure  radiation.^  The  latter 
was  done  in  reference  14  for  this  type  of  program,  and  the  former  is  presented 
herein  by  figures  7  and  8.  Also,  comparison  is  made  in  figure  8(b)  with  an  in¬ 
dependent  total  solution.  These  comparisons  give  confidence  that  convergence 
to  the  correct  solutions  is  being  obtained. 

Confidence  that  fully  converged  solutions  are  presented  is  gained  by  in¬ 
creasing  the  number  of  volume  increments  in  the  nozzle  and/or  increasing  the 
number  of  Monte  Carlo  energy  bundles  followed.  If  the  answers  do  not  change, 
which  they  did  not  for  the  cases  checked,  then  it  is  probable  that  full  conver¬ 
gence  is  present. 

It  was  found  that  20  axial  and  5  radial  increments  gave  sufficient  accu¬ 
racy.  Convergence  under  any  set  of  conditions  was  assumed  present  if  no  tem¬ 
perature  in  the  entire  distribution  changed  by  more  than  0.1  percent  between 
iterations.  This  took  three  to  four  iterations  for  most  cases. 
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start 


Assume  temperature  distribution 


(a)  Schematic  outline. 


Calculate  reference 
I  temperatures  and  properties: 


Enter  here  to  determine  radiative 
transfer  portion  of  total  nozzle 
heat  transfer 


Find  number  of  Monte  Carlo  bundles  to  account  for  sources 
due  to  flow  and  convection: 


'st,  i-l  ~  'st,  i+1, 
C„ 


Wj 

^ad‘Hw)l¥ 


If  unjxo. 

set  S(i,j)»  U(IJ), 


Go  to  beginning 
portion  of  program 


_ ^ 

Arc  all  flow  aourcc 

m 

bundles  emitted? 

Is  p  =  0? 

if  U(i,  j)  <0, 

set  S(i,  j)  ■  S(l,  j)  -  Unj) 


Did  bundle  pass 
nozzle  exft? 

Is  z>z„,v? 


Was  bundle  outside 
nozzle  as  It  passed 
exit? 

h  “  ^max 

Go  to  C  and  return 


Pick  radius  and  angles  of 

emission  at  nozzle  entrance: 

■'o ' 

^  Ff(R) 

COST?  ' 

=  (1  -  cos  27rR2 

sin  e  = 

■yiq/sini? 

Compute  Initial  Increments  and 
set  Initial  values: 
n  -  n  +  1;  Zq  -  0;  i-l; 


Since  all  flow  source 
bundles  are  now  emitted, 
emit  bundles  to  account 
for  radiant  energy  entering 
nozzle: 

p-0 


Did  bundle  leave 
through  nozzle  exit? 
Is  Z>Z,„,v? 


1$  bundle  moving 
toward  throat? 

Is  z<z„? 


Did  bundle  pass 
through  nozzle 
entrance? 

Is  z<0? 


Did  bundle  strike  wall? 
r^  ■  rad 

is  r2>[R*(z)]^? 

1^^ 

Did  bundle  strike  wall? 
r^  «  rad 

Is  r2>[R*(z)]S 

1  Find  increment  where  it  was 

Tally  bundle  In  this  gas 

absorbed: 

increment  and  re-emit  it  at 

'■(si*' 

Jmax 

new  angles  from  this  point: 
S(i,  j)  -  Sti.  j)  + 1 

_ 

ro-r,  Zo"Z 

cos  7J  - 1  -  2Rj  0  -  ZttR 

Did  bundle  pass 
through  nozzle 
entrance? 

Is  z  <0? 


Choose  angles  and  point  of  emission: 
cos  1?  - 1  -  2R;  0  -  ZttR 


Find  bundle  path  length  and  final 
axial  position  after  assumed 
absorption 

where  Kj^(i)  is  Kj^(T,  P) 
evaluated  at  conditions  in  U,  j) 
z  -  Zq  + 1  sin  sin  e 


Tally  it  in  this  Increment; 
^i  ■  Sj+i 


Was  bundle  outside  nozzle 
as  it  passed  throat? 
Za‘Zt 

Goto  C  and  return 
Is  r>[R‘(zt)f? 


Old  bundle  strike  wall  in 
this  increment? 

Zg  ■  (i  -  DAz 

Goto  C  and  return 
Is  r>  [R*(Za)]2? 


Did  bundle  strike  wall 
before  passing  entrance? 
Zg-O 

Go  to  C  and  return 
Is  r>[R*(0)]2? 


It  must  be  absorbed 
In  propellant 
Go  to  B 


Did  bundle  leave  through 
nozzle  entrance? 

Is  z<0? 


Is  bundle  traveling 

Did  bundle  strike  wall? 

toward  exit? 

r^  -  rad 

Is  Z>Zq? 

Try  next 
Increment: 

Za-Z,-Az 


Check  axial  increment  i 
to  determine  if  it  is 
where  bundle  struck 
nozzle  wall: 

Za-IAZ 


Start  checking  at  throat 
to  determine  where  it 
struck  nozzle  wall: 

Az 


Did  bundle  escape 
nozzle  end? 

Is  Z>Zmav? 


Goto  C  and  return 
Is  r>  rJ^? 

Was  bundle  outside 
nozzle  as  it  passed 
throat? 


Did  bundle  strike  wall? 

r2  .  rart 


X  -  -2 — 1  cot  n  +  r^ 
I  sine  I  0 

y  (Zg  -  Zg)  cot  e 

r.x2+y2 


Go  to  C  and  return 

Was  bundle  outside  nozzle  wall 

It  must  be 

Did  it  pass  through 

is  r>[R*(Za)]2? 

Did  it  strike  nozzle  wall? 

when  ft  passed  nozzle  end? 

absorbed  in 

.*ywV  fiozzle  end? 

^a  *  ^max 

Goto  C  and  return 

propellant 

Go  to  B 

is  Z>Z„g,? 

± 

Try  next  increment: 
Za  ■  Za  +Az 


I'  '■>[r‘(w]2? 


Tally  It  In  this  Increment 


(d)  Part  m. 

Figure  12.  -  Concluded.  Computer  flow  scheme. 
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APPEKDIX  C 


DIPFaSION  SOLUTION 

Deissler  (ref.  16)  has  derived  equations  for  the  temperature  distrihutions 
and  wall  heat-transfer  rates  due  to  radiation  in  infinitely  long  cylinders. 
Under  the  assumption  of  a  gray  gas  flowing  in  a  black  tube^  his  equation  (43) 
reduces  to 


(Cl) 


This  was  derived  using  a  second-order  diffusion  solution  incorporating  the  ef¬ 
fects  of  the  radiative  emissive  power  discontinuity  that  occurs  at  the  cylinder 
boundary.  The  Rosseland  mean  absorption  coefficient  defined  by 


must  be  used  in  this  equation. 


Use  of  the  discontinuity  in  emissive  power  at  a  solid  boundary  as  a  radi¬ 
ation  boundary  condition  is  questionable  in  the  present  case.  This  effect  ac¬ 
tually  occurs  only  in  a  model  in  which  there  is  no  conduction  of  energy  and 
wo\iLd  be  hard  to  realize  physically.  However,  the  sharp  gradient  appearing 
between  the  bulk  propellant  temperature  and  the  wall  temperature  across  the 
boundary  layer  at  any  point  along  the  nozzle  wall  provides  a  physical  situation 
at  least  approximating  this  condition,  and  the  assumption  of  a  jump  radiation 
boundary  condition  may  be  justified  here. 

This  equation  was  adapted  for  use  in  the  nozzle  by  assuming  that  the  pro¬ 
pellant  temperature  at  the  nozzle  centerline  T^,  the  absorption  coefficient 
7c,  and  the  tube  diameter  D  are  all  evaluated  at  the  axial  position  of  inter¬ 
est.  In  addition,  radiant  wall  flux  is  evaluated  on  the  basis  of  the  area  of 
each  tapering  nozzle  wall  element,  rather  than  the  cylindrical  elements  used  in 
reference  17.  For  the  equation  to  be  valid  under  these  conditions,  the  propel¬ 
lant  must  be  optically  dense  enough  that  radiation  effects  up-  and  downstream 
of  the  point  under  consideration  do  not  appreciably  affect  that  point.  This 
condition  can  only  be  tested  by  comparison  to  more  exact  solutions,  but  the 
results  will  probably  be  poor  near  both  the  exit  and  entrance  of  the  nozzle. 


The  total  heat  transfer  to  some  wall  element  i 
tive  and  convective  energy,  so  that 


hi(H  -  H^) 


will  "be  made  up  of  radia- 


(C2) 
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Loss  of  this  energy  to 
pellant  given  by 


the  wall  will  cause  a  temperature  drop  in  the  pro- 


WC  .(T.  , 


-Tp 


Substituting  ec[uations 


(C3)  and  (Cl)  into  equation  (C2)  gives 


(C3) 


Ti  = 


1-1 


(C4) 


where 


Az  JtD- 

A-  =  - - 

1  cos 


and  Az  is  the  a:xial  length  of  the  element  considered. 

Equation  (C4)  can  then  he  solved  for  if  the  inlet  propellant 

temperature,  is  given.  Tq  is  used  to  solve  for  Tg,  and  this  procedure 
is  used  to  generate  the  axial  propellant  temperature  profile.  Equation  (C2) 
is  then  used  to  determine  the  axial  wall  heat-flvix  distribution. 
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APPEKDIX  D 


HTOROGEN  PROPERTIES 

Equilibria  physical  properties  used  were  taken  from  a  number  of  sources 
and  examined  for  their  dependence  on  static  temperature  and  pressure. 

Values  of  both  enthalpy  and  radiation  absorption  coefficient  were  taken 
from  reference  17 . 

The  values  of  the  radiation  absorption  coefficient  given  in  reference  17 
are  spectral.  Examination  of  these  values  shows  a  variation  of  one  to  three 
orders  of  magnitude  over  the  spectrum  of  interest  at  any  temperature  and  pres¬ 
sure  encountered  in  the  nozzles  investigated^  with  larger  variations  in  rela¬ 
tively  low  pressiire^  low  temperature  regions.  The  Planck  mean  absorption  coef 
ficientj  defined  by 


la)  Mean  equilibrium  radiation  absorption  coefficient  of  hydrogen.  (Data  taken  from  ref.  17.) 
Figure  13.  -  Hydrogen  radiation  absorption  coefficient. 


where  e^^^  is  the  Planck  hlack  hody  energy  distrihution,  shows  a  variation  of 
many  orders  of  magnitude  over  the  range  of  either  pressure  or  temperature  ex¬ 
pected.  This  type  of  mean  results  from  consideration  of  radiation  in  gases 
that  have  negligible  self-absorption  (see  ref,  16).  The  computed  values  of  the 
Planck  mean  absorption  coefficient  are  shown  in  figure  13(a)^  and  spectral 
values  for  various  pressures  at  a  temperature  of  13^000^  R  are  shown  in  fig¬ 
ure  13(b).  Because  of  the  magnitudes  of  the  changes  with  pressure^  tempera¬ 
ture^  and  wavelength^  all  these  variables  should  be  examined.  However^  accu¬ 
rate  consideration  of  the  three  variables  simultaneously  requires  a  complex 
interpolation  routine  and  storage  of  large  amounts  of  data^  making  a  very  un- 
wieldly  program.  Variation  in  a  combination  of  two  variables  was  used^  as 
noted  in  appendix  A, 


(b)  Spectral  absorption  coefficient  for  hydrogen  gas  at 
13, 000°  R.  (Data  taken  from  ref.  17. ) 


Figure  13.  -  Concluded.  Hydrogen  radiation  absorption 
coefficient. 


Enthalpy  data  for  hydrogen  is  shown 
in  figure  14.  Since  variations  with 
pressure  are  small  in  comparison  with 
temperature  variations^  the  data  was 
approximated  by  the  dotted  lines  in  fig¬ 
ure  14.  These  are  given  by  the  rela¬ 
tions 


H  = 


3.45  T„X10"^  +  8.06 
e  s 

for  Tg  <  10,000° 


Btu/lb 
R  (Dl) 


H  ^  g  6.37  TgXlO-S  +  11.044 

for  >  10,000° 


Btu/lb 
R  (D2) 


Heat-capacity  data  are  taken  from 
reference  20  and  are  shown  in  figure  15. 
Variations  with  "both  temperature  and 
pressure  are  large,  and  neither  is  ne¬ 
glected.  Hydrogen  viscosity  data  are 
also  from  reference  20  and  are  plotted 
in  figure  16.  Variations  with  pressure 
are  small,  especially  at  high  tempera¬ 
tures,  and  the  viscosity  is  approximated 

p  =  1.293  T  X10“®  +  0.0263 
s 

lb/(ft)(hr)  (D3) 

Prandtl  number  data  from  refer¬ 
ence  20,  show  only  small  variations  with 
temperature  and  pressure  and  so  the 
Rrandtl  number  was  assumed  constant  at  a 
value  of  Pr  =  0. 65. 
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Figure  15.  -  Equilibrium  heat  capacity  of  hydrogen  gas.  (Data  taken  from  ref.  20.) 


Molecular  weight  and  ratio  of  specific  heat  data  are  found  in  references 
20  and  21.  The  molecular  weight  does  not  vary  significantly  under  the  con¬ 
ditions  expected  and,  since  it  appears  as  a  ratio  to  a  fractional  exponent  in 
equation  (A2)  its  variation  is  assumed  to  have  negligible  influence. 

The  specific-heat  ratio,  needed  to  calculate  static  conditions,  varies 
widely  with  both  temperature  and  pressure,  and  the  curves  shown  in  figure  17 
are  used. 
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